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The cosmic 21-cm background frequency spectrum related to the spin-flip transition of neutral 
Hydrogen present during and before the era of reionization is rich in features associated with phys- 
ical processes that govern transitions between the two spin states. The intervening electrons in 
foreground galaxy clusters inversely Compton scatter the 21-cm background spectrum and modify 
it just as the cosmic microwave background (CMB) spectrum is modified by inverse-Compton scat- 
tering. Towards typical galaxy clusters at low redshifts, the resulting modification is a few tenths 
milli-Kelvin correction to the few tens milli-Kelvin temperature of 21-cm signal relative to that 
of the cosmic microwave background black body spectrum. The modifications are mostly associ- 
ated with sharp changes in the cosmic 21-cm background spectrum such as due to the onset of a 
Lyman-Q radiation field or heating of neutral gas. Though low frequency radio interferometers that 
are now planned for 21-cm anisotropy measurements are insensitive to the mean 21-cm spectrum, 
differential observations of galaxy clusters with these interferometers can be utilized to indirectly 
establish global features in the 21-cm frequency spectrum. We discuss the feasibility to detect the 
spectrum modified by clusters and find that for upcoming interferometers, while a detection towards 
an individual cluster is challenging, one can average signals over a number of clusters, selected based 
on the strength of the Sunyave-Zel'dovich effect at high radio frequencies involving CMB scattering 
alone, to establish the mean 21-cm spectrum. 

PACS numbers: 98.70.Vc,98.65.Dx,95.85.Sz,98.80.Cq,98.80.Es 



Introduction — The cosmic 21-cni background involving 
spin-flip line emission or absorption of neutral Hydro- 
gen contains unique signatures on how the neutral gas 
evolved since last scattering at a redshift of 1100 and the 
subsequent reionization due to UV emission from stars 
and quasars at redshifts between 6 and 20 . The global 
signatures that are imprinted on the 21-cm background 
as a function of redshift are measurable today in terms 
of the brightness temperature frequency spectrum rel- 
ative to the black-body cosmic microwave background 
(CMB). Subsequent to recombination, the temperature 
of neutral gas is coupled to that of the CMB. At redshifts 
below ~ 200 the gas cools adiabatically, its temperature 
drops below that of the CMB, and neutral Hydrogen reso- 
nantly absorbs CMB flux through the spin-flip transition 
15,"^. At much lower redshifts, gas temperature is ex- 
pected to heat up again as luminous sources turn on and 
their UV and soft X-ray photons reionize and heat the 
gas Q. A signature is also expected from the Lyman-a 
radiation fleld produced by first sources |^ through the 
Wouthuysen-Field effect While challenging, an ob- 
servational measurement of the mean 21-cm frequency 
spectrum, in the form of modifications to the black body 
CMB spectrum, provides unique insights into the physics 
of cosmic neutral gas and the reionization history of the 
Universe 

The challenges for making a precise measurement of 
the spectrum come from the fact that it involves a 
total intensity measurement on the sky rather than 
differential measurements that are easily pursued to 
study anisotropics in the intensity. An accurate mea- 
surement of the brightness temperature spectrum re- 



quires a precise calibration of the instrument, an ac- 
counting of galactic foregrounds and point sources, 
and a control of systematics at the level of a few 
to a few tens mK 7]. Low frequency radio inter- 
ferometers that are now pursued for 21-cm observa- 
tions, such as the Low Frequency Array (LOFAR; 
http://www.lofar.org), Miluera Wide-field Array (MWA; 
http://www.haystack.mit.edu/arrays/MWA), and the 
Primeval Structure Telescope (PAST; :9j), are not sen- 
sitive to the mean spectrum as anisotropy observations 
inherently remove the mean signal on the sky. These in- 
terferometers directly measure spatial inhomogeneities of 
the 21-cm intensity, which are induced by fiuctuations in 
the neutral Hydrogen density, temperature, and veloc- 
ity, as well as spatial variations in the coupling between 
spin temperature and gas temperature such as due to 
inhomogeneities in the Lyman-a radiation field. Given 
anisotropy observations, either in the form of the three- 
dimensional power spectrum of 21-cm signal or the two 
dimensional angular power spectrum when averaged over 
a redshift interval or a corresponding frequency bin, de- 
tailed models of inhomogeneities are required to establish 
physical properties of reionization as well as the mean 
21-cm signal corresponding to that redshift or frequency 
bin |lfl|. There are large degeneracies between parame- 
ters that govern the mean brightness spectrum and those 
that lead to inhomogeneities . If signatures associated 
with the mean brightness temperature spectrum can be 
established independent of spatial variations, then some 
of the degeneracies are broken. More importantly, how- 
ever, the mean frequency spectrum alone can reveal when 
the Universe ionized and on average when the Lyman-a 



radiation field became an important source of ionization 
and wlien, if any, heating that happened to the gas. 

In this paper, we discuss the possibility to use the mod- 
ification to the cosmic 21-cm background frequency spec- 
trum by scattering via intervening electrons in galaxy 
clusters to indirectly establish global features in the mean 
21-cm spectrum generated during and prior to reioniza- 
tion. The proposed observations can be carried out with 
interferometers since the modification associated with 
low-redshift scattering can be established from differen- 
tial observations towards and away from galaxy clusters. 
Unlike an experiment to directly establish the cosmic 21- 
cm frequency spectrum at low radio frequencies involving 
a total intensity measurement on the sky, the differential 
observations with an interferometer are less affected by 
issues such as the exact calibration of the observed in- 
tensity using an external source and the confusion from 
galactic foregrounds that are uniform over angular scales 
larger than a typical cluster. The latter is the general 
case with the Galactic synchrotron background that is 
the signal at these low frequencies. The resulting modi- 
fication to the 21-cm spectrum towards a typical galaxy 
cluster at low redshifts is at the level of a few tenths mK, 
while the original spectrum has 21-cm related signatures 
with amplitudes at the level of a few tens mK relative 
to the blackbody CMB. Such a small modification chal- 
lenges an easy detection, but for upcoming interferome- 
ters, a detection could be achieved by averaging signals 
towards a sample of galaxy clusters. Such multi-cluster 
observations are easily facilitated by the fact that the in- 
stantaneous field-of-view of upcoming interferometers is 
expected to be more than 100 square degrees and one 
expects in the order of hundred or more massive clusters 
in such fields. We begin our discussion with a summary 
of photon scattering via electrons in a cluster. 
Scattering in Clusters — The inverse- Compton scattering 
of CMB photons via thermal electrons in galaxy clusters, 
generally referred to as the Sunyaev-Zel'dovich effect (SZ; 
|12l|) , is a p owerful probe of cluster physics and cosmol- 
ogy ll^lllJI' The changes to the CMB spectrum resulting 
from scattering can be studied based on the Kompaneets 
equation based on the blackbody shape of the input spec- 
trum |l5l |. The Sunyaev-Zel'dovich spectrum is that of 
a decrement and an increment, relative to CMB, with a 
cross-over at a frequency of 217 GHz; This difference in 
the spectrum relative to thermal CMB allows a variety 
of studies in multifrequency CMB anisotropy maps 0| . 

With modifications imposed to the CMB spectrum by 
neutral Hydrogen, an exact numerical calculation can be 
made to establish the low frequency spectrum emerging 
from a cluster after scattering, given an input spec- 
trum Jo(i^) = 2kBV^Tb/ in the Rayleigh- Jeans (RJ) 
end of the frequency range. Since the input spectrum 
is no longer a black-body, the standard approach for the 
CMB SZ description involving the Kompaneets equation, 
which leads to a frequency dependent term scales by the 
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FIG. 1: The scattering kernel Pi{v, vq) with vq = 100 MHz for 
a cluster with /cbTe = 7 KeV; For a different vq, kernel can be 
shifted in \og{v/vo). Any sharp features in the pre-scattered 
spectrum with variations over a frequency range less than 
the typical width of the scattering kernel will be broadened 
after scattering. The difference leads to a decrement followed 
by an increment if the incoming radiation has a dip on the 
black-body CMB spectrum, while the opposite happens if the 
pre-scattered spectrum has an increase on top of the CMB 
spectrum. 



Compton-y parameter, cannot be used here. 

Using the inverse-Compton scattering kernel, 
j/q), giving the probability that scattering has 
moved a photon from a frequency to u when the 
optical depth to scattering is r, the scattered spectrum 
can be written as (see, Ref. [3 for details) 



(1) 



The scattering probability can be calculated following 
Ref. III the case of scattering via a single electron 

moving with speed /3c, this probability is 
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where 7 = (1 - Mi = ' /^m) - '^o)]/f3iyo, the 

lower limit of the integral n\ is either —1 when 1/ < j/q or 
(j/ — i/o{l + 0))/ I3v otherwise, and the upper limit of the 
integral /j,u is either 1 when v > i/q or (:/ — fo(l — (3))/ I3v 
otherwise. 

For scattering via a population of electrons, as in the 
case of galaxy clusters, the scattering probability in equa- 
tion (2) must be averaged over the probability distribu- 
tion of electron velocities, Pe(/3): 



(3) 
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where /Sum is \v — Vfj\/{v + i'Q). We assume a thermal elec- 
tron population in clusters with an electron temperature 
Te such that the velocity distribution is given by 
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where Kn{x) is the modified Bessel function of the sec- 
ond kind. The scattering kernel P{v^ vq) is illustrated in 
Figure 1 . While we have considered the single scattering 
case, to allow for multiple scattering given the optical 
depth. If T << 1, we can write the scattering probability 
as P-r{v, uq) = (1 - t)5{v - vo) + tP{v, vq) [H m. In 
the case of CMB, the numerical evaluation of Eqs. (1) to 
(4) has been extended to study both relativistic correc- 
tions £3 and the variation associated with scattering by 
non-thermal electrons in clusters [lOl ; For the present dis- 
cussion, we ignore such complications as these only lead 
to minor corrections on top of the standard calculation. 
Results — In Figure 2, we summarize the input and the 
output spectra of the 21-cm radiation towards a cluster 
in terms of the brightness temperature that measures the 
difference with respect to CMB. Current descriptions of 
the 21-cm radiation suggest three potentially interesting 
signatures. These involve a broad dip as neutral Hydro- 
gen resonantly absorbs CMB intensity through spin-flip 
transition at a redshift range between 200 and 30 0, a 
feature related to the Lyman-a radiation field produced 
by first sources at redshifts between 30 and 20 0, and 
during reionization (at redshifts between 6 and 20), an 
increment in the brightness temperature as soft X-ray 
photons from stars and quasars begin to heat the gas |j| . 
Since the 21-cm background is formed by a line emission 
or an absorption, features that are imposed as a func- 
tion of redshift map to a unique frequency in the bright- 
ness temperature spectrum observable with telescopes to- 
day. The brightness temperature variations produced by 
above modifications to CMB by neutral Hydrogen are 
generally at the level of 20 mK to 40 mK. 

Inhomogeneities in the neutral Hydrogen distribution, 
its temperature, and the intensity of the Lyman-a pho- 
ton radiation field that couples spin temperature to the 
gas temperature via Wouthuysen-Field effect generate 
anisotropics in the brightness temperature during and 
prior to reionization. These fluctuations are at the level 
of a few mK to tens of mK 0| . A measurement of these 
anisotropics, in the form of an angular power spectrum 
is now considered to be one of the main goals of planned 
low-frequency interferometers. Additional fluctuations 
are also generated during the propagation such as due to 
gravitational lensing and as we discuss later through 
scattering in galaxy clusters. The intrinsic anisotropics, 
however, are dominated by those related to low radio fre- 
quency foregrounds and techniques are now developed to 
remove confusing sources |22j . 

The brightness temperature difference towards and 
away from a cluster is shown in Figure 2(b). First, note 
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FIG. 2: (a) The incoming radiation spectrum plotted in terms 
of the brightness temperature relative to CMB. Three features 
are shown: A broad absorption due to collisional coupling be- 
tween spin and gas temperatures when 30 < z < 20oQ, cou- 
pling induced by Ly-a radiation from first luminous sources 
when 20 < 2 < 30 p (labeled Ly-a), and an increase in the 
gas temperature when heated by soft X-rays during reion- 
ization at a redshift below 20 ,4j (labeled X-ray), (b) The 
difference spectrum after scattering through a galaxy cluster 
with an electron temperature of 7 keV and an optical depth 
of 0.005, consistent with known clusters at low redshifts |2^ 
with the optical depth profile smeared over a beam of 5 ar- 
cminutes to be consistent with typical beam sizes expected for 
upcoming interferometers. The overall offset in the spectrum, 
at a brightness temperature ~ -0.35 mK with no frequency 
dependence at the Rayleigh- Jeans tail, is due to scattering of 
the CMB black-body with no 21-cm signatures. The 21-cm 
signatures lead to variations up to 0.2 mK, While the broad 
dip at very low frequencies leads to a decrement and an in- 
crement, Ly-Q and X-ray features combine to given an overall 
increase towards clusters at a frequency around ~ 75 MHz; 
This is associated with the sharp dip related to the Lyman-a 
induced coupling of gas and spin temperatures at a redshift 
around 20. The dot-dashed line is the scattered spectrum only 
involving the broad dip and the X-ray heating increment; the 
latter leads to an increment followed by a decrement. 



the overall amplitude difference between the two pan- 
els. While the brightness temperature in the unscattered 
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FIG. 3: The brightness temperature difference spectrum to- 
wards and away from a cluster with a beam-smeared aver- 
age optical depth of 0.005 and an electron temperature of 7 
keV that can be targeted with upcoming low-frequency radio 
interferometers (same as Figure 2b); the optical depth and 
electron temperature values are typical of low redshift clus- 
ters |23| . Since underlying 21-cm signatures are global, one 
can also stack spectra towards a large sample of clusters to 
improve the signal-to-noise ratio. The errors bars assume ex- 
perimental parameters similar to SKA with a system noise 
temperature of 1000 K at 100 MHz, and a scaling of to 
other frequencies, a 5 MHz bandwidth, synthesized beam of 
5', and a year-long integration. The larger errors indicated 
by horizontal lines assume observations towards an individ- 
ual cluster, while error bars show the improvement by av- 
eraging over observations of 10 clusters, though in planned 
observations ^ few hundred massive clusters will be imaged 
in a given field-of-view instantaneously. A large uncertainty 
in establishing the exact signal-to-noise for detection is the 
foreground distribution within clusters, such as point sources, 
and the extent to which such sources can be removed |23I. 



spectrum is generally at a level of 30 mK, the bright- 
ness temperature spectrum after scattering is an overall 
shift of about -0.35 (t/0.005) mK, independent of the 
frequency, with variations of order 0.05 mK to 0.2 mK 
imprinted on it. As an example here, we selected a cluster 
with an optical depth r of 0.005 and an electron temper- 
ature of 7 keV. These values are consistent with known 
electron distribution of low-redshift clusters {23] when 
the optical depth profile is averaged over a beam size of 
5 arcminutes. The 5 arcminutes beam size is the typi- 
cal of what is expected for upcoming low frequency radio 
interferometers at a frequency of 100 MHz |22| . 

The overall offset in the mean temperature to a value of 
—0.35 mK is the usual SZ effect associated with the black- 
body CMB spectrum at the Rayleigh- Jeans part of the 
frequency spectrum. On top of this signal one finds 21- 
cm signatures in the form of decrements and increments. 
The broad absorption in the redshift range between 200 
and 30 modifies to a sharp decrement of 0.05(r/0.005) 
mK followed by an increment of < 0.05(t/0.005) mK. 



While such very low frequencies will not be observable 
with upcoming radio interferometers, in the frequency 
range easily accessible to first generation 21-cm interfer- 
ometers is the signature related to reionization. Since 
the gas heating increment quickly follows the Ly-a re- 
lated dip, the overall change after scattering in clusters 
is a decrement followed by a sharp increment and another 
decrement. The latter signatures are model dependent as 
they require assumptions on when sources turn on, the 
intensity of the Lyman-a radiation produced by these 
sources, and the intensity of the soft X-ray background 
that heats the gas. We also considered the case where Ly- 
a signature is not present, but only a signature related to 
gas heating (dot-dashed line in Figure 2b). In this case 
modifications to the spectrum have smaller amplitudes, 
but is an increment followed by a decrement. Clearly, the 
amplitudes and the cross-over frequencies are sensitive to 
exact details of the reionization process and a measure- 
ment of the spectrum in Figure 2(b) should be a goal for 
upcoming experiments. 

If the pre-scattering brightness temperature spectrum 
has signatures at the level of 30 mK to 40 mK over a 
narrow range in frequency, as in the Ly-a signature of 
Figure 2(a), resulting signatures in the post-scattering 
brightness temperature spectrum towards clusters are 
potentially measurable with current and upcoming in- 
terferometric arrays. To illustrate this possibility, in Fig- 
ure 3, we plot the expected errors on the brightness tem- 
perature assuming an array with a system temperature. 



Tsys, of 1000 K at 100 MHz, with T^. 



a band- 



width Aj/ of 10 MHz for observations, an integration time 
of 1 year, and a collecting area that of the Square Kilo- 
meter array (SKA) with A ^ 1km?. For interferometric 
observations, the pixel noise in the synthesized beam is 
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Since the background spectrum is the quantity of inter- 
est, given prior knowledge on the optical depth and elec- 
tron temperature, one can appropriately scale and com- 
bine observations towards multiple clusters to increase 
the signal-to-noise ratio. In this case, we assume the im- 
provement in AT with \fW^\ of the number of clusters 
A^ci used. 

The errors shown in Figure 3 are those using a single 
cluster and the improvement by averaging over signals 
towards 10 clusters. When integrated over the whole 
spectrum, relative to the case where the scattering effect 
is due to CMB alone, one detects variations at the level 
of 2 a towards a single cluster, but better than 10 a 
when averaged over ^ 10 clusters. For smaller arrays 
that will be soon be operational, such as MWA with an 
area 10 times smaller than the SKA, due to the large in- 



5 



stantaneous fields-of-view of 300 deg.^ [2j|, ~ 10^ clusters 
will be imaged and a better than 5 a detection of vari- 
ations related to Ly-a coupling or X-ray heating can be 
obtained by stacking signals towards ~ 100 clusters. The 
stacking procedure requires prior knowledge on the clus- 
ter optical depth profile and the electron temperature. 
This can be obtained a priori based on SZ observations 
related to the CMB spectrum alone at high radio fre- 
quencies. In the near future, all-sky maps from Planck 
will provide a substantial catalog of massive SZ clusters 
with profiles averaged over the Planck's angular scale of 
5 arcminutes Since the angular scale naturally cor- 
responds to those of low-frequency observations, one can 
use information on the high frequency SZ profile as a 
way to normalize low frequency signals when stacking to 
establish signatures related to the 21-cm spectrum. 

The calculation related to the modification resulting 
from scattering in clusters assume that the incoming 
spectrum through the cluster is same as the 21-cm fre- 
quency spectrum outside the cluster. Any difference be- 
tween the two will lead to an additional signal in the dif- 
ference towards and away from the cluster. While such an 
extra difference can aid the detection, it can also com- 
plicate an easy interpretation of the measurements es- 
pecially when attempting to establish global signatures 
in the mean 21-cm spectrum as one must account for 
arcminute-scale anisotropy in the incoming signal itself. 
By carefully averaging signals over a large number of clus- 
ters, spatial variations in the incoming signal may be ac- 
counted for. This scattering related differences towards 
and away from clusters may also be used as a probe of 
the typical size scale of reionization, such as due to in- 
homogeneities in the Lyman-a radiation field from first 
UV sources. 

Studies towards a large sample of clusters is possible 
since upcoming low-frequency interferometers will have 
instantaneous fields of view of order 300 deg.^ or more. 
One would naturally expect of order a few hundred clus- 
ters in such fields with adequate mass to produce measur- 
able differences. The same clusters will also be detectable 
at high radio frequencies with CMB experiments such as 
Planck. The latter observations are useful to find suitable 
cluster candidates for low-frequency studies based on the 
strength of the SZ signal at high frequencies and an es- 
timate of the optical depth profile, in combination with 
X-ray estimates of the temperature. The SZ spectrum 
alone, at high-frequencies, may provide independent es- 
timates of the temperature (at the level of a few keV), 
which may be useful for low-frequency cluster studies in 
the absence of X-ray data 0| . 

Before concluding, we note that for all these stud- 
ies, the main worry is complications resulting from fore- 
grounds. Since observations are differential any smooth 
foreground at angular scales more than the cluster will 
be removed. This will remove a significant fraction 
of the Galactic foreground signal. The main worry is 



point sources within clusters and a potential smooth syn- 
chrotron background associated with electrons in clusters 
that appear as radio halos and radio relics . The fore- 
grounds, however, are expected to be controlled through 
a combination of direct point-source removal, when re- 
solved, and through spectral information; Techniques 
that are developed to remove foregrounds in anisot ropy 
data can be easily applied for cluster observations |22l |. 
Unfortunately, estimates on the level of confusion within 
clusters are highly uncertain due to the lack of cluster 
observations at low radio frequencies. Upcoming arrays 
will provide a wealth of data to further understand the 
expected level of confusion. 

To summarize, the cosmic 21-cm background frequency 
spectrum related to the spin-flip transition of neutral 
Hydrogen present during and before the era of reioniza- 
tion is rich in features associated with physical processes 
that govern transitions between the two spin states. The 
intervening electrons in foreground galaxy clusters in- 
versely Compton scatter the 21-cm background spectrum 
and modify it just as the cosmic microwave background 
(CMB) spectrum is modified by inverse-Compton scat- 
tering. While low-frequency interferometer arrays are 
built or planned to primarily study 21-cm anisotropics, 
they can easily target galaxy clusters to establish signa- 
tures in the 21-cm background. We strongly encourage 
an analysis of signals towards and away from clusters as 
a way to establish signatures related to Lyman-a radi- 
ation intensity field or gas heating during and prior to 
reionization. The analysis is challenging as the expected 
signal difference due to scattering is at most 0.2 mK for 
low redshift clusters with typical optical depths of order 
0.005 and may be further complicated by the underlying 
anisotropics associated with topology of reionization it- 
self. The availability of imaging data for a large numbers 
of clusters (~ 10^) instantaneously with low frequency ra- 
dio interferometers for 21-cm studies as well as the avail- 
ability of SZ data for such clusters soon with experiments 
like Planck will allow joint studies to be conducted. 
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